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We compute the order o? a corrections to the one particle inclusive electroproduction cross section 
of hadrons with non vanishing transverse momentum. We perform the full calculation analytically, 
and obtain the expression of the factorized (finite) cross section at this order. We compare our 
results with HI data on forward production of 7r°, and discuss the phenomenological implications of 
the rather large higher order contributions obtained in that case. Specifically, we analyze the cross 
section sensitivity to the factorization and renormalization scales, and to the input fragmentation 
functions, over the kinematical region covered by data. We conclude that the data is well described 
by the 0(ct 2 a ) predictions within the theoretical uncertainties and without the inclusion of any physics 
j /^ ' content beyond the DGLAP approach. 
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Introduction 



7-H ! The precise measurement of final state hadrons in lepton nucleon deep inelastic scattering constitutes an excellent 
i benchmark for different features of perturbative quantum chromodynamics (pQCD). These processes are crucially 
sensitive to the three main ingredients of pQCD: the parton content of the nucleon, the hadronization mechanism 
■ of partons into the detected final state hadron, and the parton radiation before and after the interaction with the 
O |' electromagnetic prove. 

The first of these ingredients is well characterized by modern parton distribution functions (PDFs). The knowledge 
on these distributions have become increasingly precise as a result of two decades of high precision inclusive mea- 
surements, and the corresponding QCD analyses, driven by the role of PDFs as inputs for theoretical predictions for 
any experiment involving initial state hadrons 0,13 ■ Although the high degree of accuracy attained by PDFs, less 
. £^ \ inclusive observables, sensitive to flavor combinations of PDFs other than those relevant in inclusive measurements, 
improve the insight and provide a further check on the universality of PDFs and on factorization. 

The second ingredient is addressed by the so called fragmentation functions (FF), which are rapidly evolving 
following the path of PDFs, but without attaining yet the refinement of the latter |3, 0. Most of the data used 
to determine these distributions, which come essentially from electron-positron annihilation into hadrons, give no 
information on how the individual quark flavor fragments into hadrons and leave a considerable uncertainty on the 
gluon density. For these reasons, the NLO analysis of one particle inclusive data is crucial for the extraction of 
fragmentation functions. 

The last ingredient concerns higher order QCD calculations, which have been explored and validated for most 
processes up to next to leading order (NLO) accuracy, and are currently being extended even beyond that point. 
For the one particle inclusive processes only very recently there has been progress beyond the leading order (LO) 
jHElElElil- However, up to now there were no analytic computation of the 0{a 2 s ) corrections for the electroproduction 
of hadrons with non vanishing transverse momentum. The analytic computation of the 0(pt 2 ) corrections allows us 
to check factorization in a direct way, which means that collinear singularities showing up in the partonic cross 
section factorize into PDFs as required by inclusive deep inelastic scattering, and into FFs for electron-positron 
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annihilation into hadrons. As a consequence of this explicit cancellation, the resulting cross section is hnite and can 
be straightforwardly convoluted with PDFs and FFs in a faster and more stable numerical code, compared to what can 
be usually obtained in numerical implementations using either the subtraction Q or the slicing 0, 0] methods. The 
analytical result is still sufficiently exclusive and keeps the dependence on the rapidity and the transverse momentum 
of the produced hadron, allowing a detailed comparison with the experimental data. 

In this paper we compute the order a 2 corrections to the one particle inclusive leptoproduction cross section of 
hadrons with non vanishing transverse momentum. We perform the cancellation of collinear singularities analytically, 
and obtain the full expression of the factorized, and thus finite, cross section at this order. The outline of the paper 
is as follows: in the next section we summarize the relevant kinematics and details about the phase space integration 
for the 0(a 2 ) contributions to the cross section, together with the conventions and notation adopted. In section II 
we compute the corresponding real and virtual amplitudes, we discuss the nature of the singularities that contribute 
to them, we analyse the factorization of collinear singularities, and pay special attention to the scale dependence 
induced in the cross section by this factorization procedure. Section III deals with the phenomcnological implications 
of the new corrections. Specifically, we compare our results with data on forward production of 7r°, presented recently 
by the HI collaboration, as an example, and evaluate the phenomenological implications of the rather large higher 
order corrections. Special attention is paid to the cross section sensitivity to the factorization scale choosen, and to 
the fragmentation functions input as sources of theoretical uncertainties. We also analyze the consequences of the 
forward pion selection on the LO and NLO underlying partonic processes, finding this kinematical suppression as the 
main reason for the unusual difference between the LO and NLO estimates. In agreement with the results obtained 
in 0, we conclude that the data is well described by the pure Dokshitzer-Grivov-Lipatov-Altarelli-Parisi (DGLAP) 
0{a s ) predictions within the theoretical uncertainties, but without need to appeal to resolved photon contributions, 
as suggested in Q. 



I. KINEMATICS 



We begin with the kinematical characterization of the one particle inclusive deep inelastic scattering processes. 
Since the choice of variables re quir ed to deal with the singularity structure of electroproduction is different from those 
used in both photo-production |10L IllL and electroproduction in the very forward region 0, 0] , in the following we 
discuss it in some detail. We consider the process 



l(l) + P{P) — > Z'(Z') + h{P h ) + X, 



(1) 



where a lepton of momentum I scatters off a nucleon of momentum P with a lepton of momentum I' and a hadron h 
of momentum Ph tagged in the final state. Omitting target fragmentation at zero transverse momentum, which has 
been discussed at length in 0, [(J , the cross section for this process can be written as 



dx b dQ' 



MOD h/j (() 



da 



(n) 



dx B dQ 2 dPS (n) 



(2) 



where <jf^ is the partonic level cross section corresponding to the process 

1(1) + i(pi) — > + + n_1 additional partons , 



(3) 



before renormalization of the coupling constant and factorization of collinear singularities. and D^/^CC) are the 

bare parton densities and fragmentation functions, and dPS^"-* the n-parton phase space. £ is the proton momentum 
fraction carried by the parton i and £ is the fraction of parton j momentum taken away by the final state hadron. In 
addition to the usual DIS variables, 



Q 



-{l'-lf 



x B 



Ve = 



2P-q 

we define Mandelstam variables both at parton and at hadron level: 

s={q+ Pi ) 2 
t = —2 q ■ kj 
u = —2 pi ■ kj 



P_q 
P-l ' 



S H = (P + l) 2 , 



S=(q + P) 2 , 
T=-2q-P h , 
U=-2P-P h . 



(4) 



(5) 
(6) 
(7) 
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respectively. The above definitions imply 
fl = ££-Q 2 (l-0, 



T 



u=-U. 

c 



(8) 



Notice that the t > (T > 0) region exists only for Q 2 ^ 0, feature that considerably reduces the integration region 
in the case of photo-production. The following step is the definition of suitable partonic variables to characterize 
the phase space. The choice of these variables is critical for the identification and further prescription of collinear 
singularities in the partonic cross section. We find particularly useful the variables 



y 



(Q 2 + s)(s + t + u) 



s s{Q 2 + s + u) ' 

with y, z € [0, 1]. In terms of these partonic variables, the n-particle phase space can be factorized as: 



dPS (n) =dPS ( " 1] dydz, 



(9) 



(10) 



>-i) 



where dPS includes the phase space of the 'spectator' partons (those that not fragment into the detected final 
state hadron) and the corresponding jacobian. For example for n = 3, in D = 4 + e dimensions we have 



dPS (3) = (*) 
V47T/ 



(4 7r) 4 r(l + e) 
x sin 1+e /?i sin £ f3 2 dfli dfi 



(l-y)z^ 2 f/ 2 (l-yY(l-zy/ 2 dydz 



(11) 



where f3\ and (3% are the angles defined by the spectator partons in their center of mass frame. In terms of the 
factorized phase space, eq.@ reads 



dx b dQ 2 



E 



i j-i pi i-i 

d£ / dC, \ dy I dz 
i) Jo Jo Jo 



da 



(n) 



dxB dQ 2 dy dz 



(12) 



where da^jdx^ dQ 2 dy dz is the partonic cross section already integrated over the spectator partons, and with the 
adequate normalization. Finally, changing variables from (£,£) to hadronic transverse momentum pt and rapidity 77, 
defined in the center of mass frame of the proton and the virtual photon, we find 



da h 



E 



iin \ PT \(Q 2 + S) J eV \j^ l-yJo 



c 2l 

l+e^i dy 



dxB dQ 2 dp^ dr\ 

In terms of the hadronic variables, ^ and C are given by 

Q 2 (l-y)(l-z) + Sye- 2 ri 



dz 



da 



(n) 



dxB dQ 2 dy dz 



(Q 2 +S)(1- y) (1 - z) 



C 



e^pi 



(13) 



(14) 



Clearly, the transformation is singular at y = 0, y = 1 and z = 1, however these points are excluded by \pr\ > 
(notice that \r)\ is bounded from above), as can be seen from the limits in eq. i|13|) . Finally, in order to obtain more 
compact expressions for the partonic cross sections, it turns out to be convenient to introduce the auxiliary variable 



Xb 



Q 2 



£ Q 2 



(15) 



II. ORDER a a AND ai PARTONIC CROSS SECTIONS 



The partonic cross sections in ea. (|13f) are calculated order by order in perturbation theory and are related to the 



parton-photon squared matrix elements H for the i + 7 — » j + X processes 



da 



(n) 



dxB dQ 2 dy dz 



£x B S H 



4x 



(16) 
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Matrix elements are averaged over initial state polarizations, summed over final state polarizations, and integrated over 

- — 

the spectator partons (i.e. integrated over dPS ). a em stands for the fine structure constant and e is the electron 
charge. Finally, Ym and Yl are the standard kinematic factors for the contributions of each photon polarization and 
are given by, 



Y, 



M 



2y 2 



Y L 



l + A{l-y e ) + {l-y e f 



(17) 



The first contribution to the cross section comes from the partonic tensor at order a s , as, in the naive parton 
model (©(a^)), final state hadrons can only be produced with \px\ — in the proton-virtual photon rest frame. At 
order a s , the partonic cross sections have no collinear divergences provided \pr\ > 0. Up to order e, they are given by 



da, 



(i) 
</</ 



c q C e 



dx b dQ 2 dy dz £ x B Sjj 



C F <Y, 



M 



(g + y) 2 + 2(1- g-y) (,,e\ 



(9 - v) 



2(1-0) (l-y) 



da, 



(i) 

19 



dx b dQ 2 dy dz £ x 2 B S B 



da 



(1) 



Y L 



CqCe cAy m 



dx b dQ 2 dy dz £ x 2 B Sjj 



CqCe tAy 



M 



2gy(l + -L 1 



+ 0{e 2 )\5{z), 



+ Y L 



Yl 



[l- e )y \ 2 V + 2 (l-g)y 
2g(l-y) + +0(e 2 )\6(z), 



(18) 



(19) 



((1-g) g+(l-y) y) 



(1-2 (1-g) g-2 (1-y) y) / e 

(1-1/) 1/ V 2 V ' ~ CI // 

4 (1-g) q (l + |ii) -2e (1-g) g] +G(e 2 )U(z) 



(20) 



where 



c, = a 2 m 27r (2 + e) e 2 



_a s ( Q 2 ^ 6/2 
Ct ~ 2^ /r (,4^ 



/r = 



r(l + e/2) 

r(i + e) 



(21) 



and 



L\ = log 



/(l- e )(l-y)y 



(22) 



At order-a 2 , the partonic cross sections receive contributions from the following reactions: 



Real contributions 



Virtual contributions 



( i + q(q) -> g + g + q(q) 
7 + q% (qi) -> qi (<ii) + Qj + qj (* ^ i) 

7 + ft (ft) qi (ft) + ft + ft 

I 7 + ff -> g + q + q 
i + q{q) -> g + q(q) 
i + g q + q 



(23) 



where any of the outgoing partons can fragment into the final state hadron h. Order a 8 and a 2 contributions in the 
very forward region and its singularity structure have already discussed in deep in references 0, ■ In this section 
we analyze the \pt\ > region, and examine the nature of the singularities that it involves. These contributions are 
computed in d = 4 + e dimensions, in the Feynman gauge, and considering all the quarks as massless. Al gebr aic 
manipulations were performed with the aid of the program Mathematica [l3| and the package Tracer [lj] to 
perform the traces over the Dirac indices. 

The order a 2 partonic cross sections can be obtained from the corresponding quark and gluon initiated amplitudes 
as in references ja, |(| , taking care of the appropriate flavor discrimination. The angular integrations can be performed 
with the standard techniques flU Il6j | , taking into account the additional complications of the one particle inclusive 
case: the necessity of collecting to all orders the potentially singular factors in the three particle final state integrals. 
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For the integrals that are known to all orders in e, this is not a problem, while for those which are only known up 
to a given order a careful treatment is required. Once the angular integrals are performed, matrix elements are still 
distributions in three variables, y and z and g, regulated by the parameter e. 

At variance with the very forward case (\pt\ — 0), where the integration over final states leads to overlapping 
singularities along various curves in the residual phase space, here the only remaining singularities are found at z = 
and thus they can be dealt with the standard method. After combining real and virtual contributions to a given 
partonic process, the cross section can be written as 



da 



(2) 



r 2 



dxg dQ 2 dy dz £ x 



B H 



v[%(g,y,z) +a 



(2) 



(g,y,z) + 0(e) 



(24) 



(2) (2) 

where the coefficient of the single poles, V\l^g,y, z), as well as the finite contributions C\j (g,y, z), include 'delta' 
and 'plus' distributions in z. The IR double poles present in the individual real and virtual contributions cancel out in 
the sum, providing the first straightforward check on the angular integration of real amplitudes and the loop integrals 
in the virtual case. In the real terms, the above mentioned double poles come from the product of a pole arising in 
the integration over the spectators phase space (i.e. integration over pi and /?2 in (f 1 If) ^ and a single pole coming from 
the expansion of z~ 1+e factors. Double poles in the virtual contributions always arise from loop integrals. 

The remaining singularities, contributing to the single pole, are of UV and collinear origin. The former are removed 
by means of coupling constant renormalization, whereas the latter have to be factorized in the redefinition of parton 
densities and fragmentation functions. The redefinition of parton densities is exactly the same as in totally inclusive 
DIS whereas fragmentation functions are renormalized as they are in one-particle inclusive electron-positron anni- 
hilation. Typical expressions for renormalized parton densities and fragmentation functions, up to order a 2 and in 
the MS factorization scheme, can be found, for example, in Refs.^?J an d EH respectively. Factorization of collinear 
singularities and cancellation of the UV ones, then impose 



V[%(g,y 7 z) = 2(C l 



(i) 



>(o) 



(25) 



where the C^} correspond to the finite (0(e )) terms in the 0(a s ) partonic cross sections of eqs 



(J2DJ. Pij are the standard LO Altarelli-Parisi kernels, and 



(PI). HU and 

denotes the appropriate convolution coming from the 



factorization recipes. 

The factorized, and thus finite, partonic cross sections have terms proportional to S(z), terms containing 'plus' 
distributions, and purely functional contributions. The logarithmic 'plus' contributions have their origin in the 
multiple emission of soft-gluons and can therefore be predicted by taking the order a s expansion of the corresponding 
resumed cross-section. For a partonic subprocess initiated by a parton i, where a parton j fragments, and with a 
gluon and a parton s as spectators, 7 + i — ► g + j + s, the result is 



da 



(2) 



dxB dQ 2 dy dz 



da 



(i) 



dxB dQ 2 dy dz 



(26) 



where the general color factor Ck corresponds to Cf if k is a quark and to Ca if it is a gluon. The agreement with 
this prediction provides a further test on our results. 

Since the factorized coefficients get contributions from both the real and virtual process at order a 2 , together with 
finite terms coming from the renormalization and factorization procedure, their explicit expressions are considerably 
long and thus are omitted here (24j. 

Notice that the renormalization and factorization processes introduce scale dependent terms in the final cross section 
which partially cancel the scale dependence induced by the coupling constant, parton densities, and fragmentation 
functions. The structure of these terms follows that of the factorization contributions in ea. (|25|l . 



C$ (8 Pli In 



Ml 
Q 2 



C 



(i) 



Pjk In 



Mf, 

Q 2 



■ In 



Ml 

Q 2 



(27) 



where Mr, is the renormalization scale and Mf and Md are factorization scale for parton densities and fragmentation 
functions, respectively. 

In order to visualize the magnitude of the higher order corrections, in Figure 1 we show the ratio between the 
order a 2 and a s cross sections for tt° production (K- factor) as a function of pr for Q 2 — 200 GeV 2 , integrated over 
rapidity and for different values of in in the kinematically allowed range. As input parton densities and fragmentation 
functions we choose the MRST02 and KKP NLO and LO sets, respectively. In the following we refer to the 
convolution of 0(a 2 ) cross sections and NLO densities as NLO prediction, whereas 0(a s ) cross sections convoluted 
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FIG. 1: K-factor 


as a function of pr- 



with LO densities define the LO estimate. The renormalization and factorization scales Mi Mj, and M|, were taken 
to be the average between the two main physical scales of the process, namely the transverse momentum of the final 
state particle and the virtuality of the photon as 



Ml 



Ml 



Q 2 



Pt 



(28) 



The K-factor exhibits the characteristic behavior of higher order corrections; they increase at low transverse momentum 
and also at low xb- At very high px, where the LO estimate becomes larger than the NLO prediction, threshold 
effects become dominant and the perturbative expansion at fixed order in the coupling constant is not expected to be 
reliable. 

The dependence on the particular choice for the factorization and renormalization scales is expected to be weaker 
at NLO than at LO. In Figure [21 we show this dependence plotting the rate between the cross section evaluated at 
an arbitrary scale /i 2 and the cross section at [Aq = (Q 2 + p^)/2 as function of the rate ^ 2 //j,q. As in the previous 
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FIG. 2: Scale dependence of the electroproduction cross section. 



plot Q 2 = 200 GeV 2 , but xb = 0.01 and we integrate over the allowed range, starting from \pt\ > 3.5 GeV. As 
expected, the scale dependence is milder for the NLO estimate, although it is not negligible. 

Notice that our NLO estimate focus on the 'direct' coupling of photons to partons, without taking into account the 
'resolved' photon contributions, as those_computed with virtual photon parton densities. These contributions have 
been carefully analyzed in references 
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III. PHENOMENOLOGY 

Recently the HI [20| collaboration has presented an improved measurement of the production of neutral pions in 
collisions between 27.6 GeV positrons and 820 GeV protons. Neutral pions are required to be produced within a small 
angle 6 n from the proton beam in the laboratory frame [9„ € [5°, 25°]), with an energy fraction x^ = E^/Ep > 0.01 
and 2.5 < pr < 15 GeV. The data confirmed previous measurements which suggested that QCD LO predictions 
underestimate the cross section at low xb HH]. On the other side, predictions based on BFKL dynamics [2^], or on a 
large virtual photon content |2.'"{l| seemed to provide better descriptions. 

The disagreement between the HI data and estimates based on 0(a s ) cross sections convoluted with LO parton 
densities and fragmentation functions can be as large as an order of magnitude, depending on the particular kinematical 
region. This discrepancy is far larger than the typical K-factor found in the previous section, what suggests the onset 
of a physical mechanism different to leading or next to leading order DGLAP dynamics. 

However, several non-negligible effects are present at the particular kinematical regime of the experiment, which 
are responsible for a large difference between the LO and NLO estimates. The first one is the stringent cut on the 
pion production angle in HI data, which suppresses LO and NLO contributions in a different way. The suppression 
of LO configurations is proportionally bigger than for NLO, implying an effective K-factor much larger than the one 
found for the cross section without cuts. The second important feature is the rather low value of the scales involved 
(pt and Q 2 ) which enhance the uncertainty due to the particular choice for the factorization scale, even in the NLO 
calculation, as it has been pointed in pfl. This is particularly significant for the lowest Q 2 bins. Finally, there is also a 
large uncertainty factor in the theoretical prediction coming from fragmentation functions. Although fragmentation 
functions reproduce fairly well e + e _ annihilation into hadrons, they show large differences when they are used to 
compute deep inelastic semi-inclusive cross sections. 

In Figures and 0] we show the LO and NLO predictions for the electroproduction of neutral pions as a function 
of xb and pt, respectively, in the kinematical range of the HI experiment, together with the most recent data for 
the range pt > 3.5 GeV. The cross sections are computed as described in the previous sections, applying HI cuts 
and using MRST02 parton densities 0] . Similar results are found using other sets of modern PDFs. For the input 
fragmentation functions, we use two different sets, the ones from reference Q denoted as KKP and those from £| 
referenced as K. We set the renormalization and factorization scales as in eq. (27) and we compute a s at NLO(LO) 
fixing Aqcd as in the MRST analysis, such that a s (M z ) = 0.1197(0.130). 
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FIG. 3: LO and NLO cross sections, including experimental cuts as explained in the text, as a function of xb- HI data (2{j for 
the range pr > 3.5 GeV are also shown. 



The plots clearly show some of the features mentioned above. On the one hand, the NLO cross sections are much 
larger than the LO ones, even by the required order of magnitude in certain kinematical regions, once the forward 
7T° selection applied by HI is implemented. The position of the maximum for the Xb distribution is also shifted to 
lower xb values, agreeing nicely with the experimental shape. Cross sections differential in p? show similar features, 
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FIG. 4: Cross section as a function of pr, data and cuts as in Figure 



however the difference between LO and NLO decreases as pr increases. 

The uncertainty due to the choice of a fragmentation functions set is also quite noticeable, this fact driven by the 
different gluon content of the two sets considered here. Low Q 2 bins seem to prefer KKP set, which have a larger 
gluon-fragmentation content, whereas for larger Q 2 both sets agree with the data within errors. LO estimates show 
a much smaller sensitivity on the choice of fragmentation functions, since gluon fragmentation does not contribute 
significantly to the cross section at this order. 




FIG. 5: NLO cross sections as a function of xb- The central (solid) line corresponds to setting the factorization and renormal- 
ization scales to (Q 2 +Pt)/2 and the upper and lower (dashed) ones to (Q 2 + Px)/4 and Q 2 + Pt respectively. 



As mentioned, the dependence of the cross section in the choice for the renormalization and factorization scale is 
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also an important source of uncertainty even at NLO. In Figure |5] we show the NLO prediction with the standard 
choice for the scale, MRST02 parton densities and KKP fragmentation functions (solid line) as in Figure |3 together 
with HI data and the estimates with a scale twice as large (lower dashes) and another scale half of the former (upper 
dashes) . 

Finally, in order to illustrate the effects of the forward selection criteria, in Figure we show the effective K-factor 
for the lowest Q 2 bin, with and without taking into account the constraints on 0^ and x w . Although the low values 
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FIG. 6: K-factors for the lowest Q 2 bin of Figure 01 as a function oi xb with and without the experimental cuts. 

of Xb and Q 2 lead to a very large K-factor, the forward selection, typically enhances it by a factor of three. Notice 
that the process 7 + g — > g + q + q becomes active at 0(a 2 ) and indeed turns out to be responsible for most of the 
correction, as it is illustrated in Figure[7| where we show the different contributions to the cross section discriminated 
by the underlying partonic process. 

The rather large size of the K-factor can, then, be understood as a consequence of the opening of a new dominant 
('leading-order') channel, and not to the 'genuine' increase in the partonic cross section that might otherwise threaten 
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FIG. 7: Contributions to the cross section discriminated by the underlying partonic process for the lowest Q 2 bin of Figure |3 
including experimental cuts. Processes qg and qq give negligible contributions and are not shown in the plot. 
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perturbative stability. The dominance of the new channel is due to the size of the gluon distribution at small xb and 
to the fact that the HI selection cuts highlight the kinematical region dominated by the 7 + <7 — > <7 + <Z + <? partonic 
process. In particular, without the experimental cuts for the final state hadrons, the gg component represents less than 
25% of the total NLO contribution at small xb, which is dominated by the gq subprocess. The forward selection is also 
responsible of the scale sensitivity of the cross section, as it supresses large components with small scale dependence 
whereas it stresses components as gg whose scale dependence would be partly canceled only at NNLO. 



IV. CONCLUSIONS 



We have presented the analytical calculation of the differential cross section for semi-inclusive production of a 
hadron, with non vanishing transverse momentum, in DIS at next-to-leading-order in QCD. As for any semi-inclusive 
process, the necessity of integrating the phase space of the unobserved particles, but keeping the full dependence on the 
variables characterizing the final state hadron (and thus of the parton from where it comes), makes the computation 
of higher order corrections much more involved than the inclusive case. In the present case we showed that, with 
a suitable parameterization of the phase space, the necessary integrations can be performed analytically and the 
remaining singularities can be dealt with standard prescription recipes, without the need of substraction or phase 
space slicing methods. 

We found that the order a 2 s corrections are important, leading to large K-factors. The main contributions to these 
corrections come from the partonic subprocess "f + g—>g + q + q which appears for the first time at that order. The 
appearance of new channels also leads to quite a significant factorization scale dependence even at the NLO level. 

Concerning the phenomcnological consequences of our results, we compared them with recent data coming from 
the HI experiment at HERA [2(|. Within the uncertainties arising from the scale dependence and the particular 
sensitivity of the results to the gluon hadronization mechanism, parameterized in the fragmentation functions, we 
found a very good agreement between data and theoretical expectations for both the xb and px distributions. In 
particular, fixing the factorization and renormalization scales as just the average between the photon virtuality and 
the transverse momentum of the final state hadron, both distributions are well described by purely DGLAP evolution. 
We also found that the experimental cuts applied to the HI data play a crucial role, boosting the NLO corrections, 
and thus explaining the unusual poor description of the LO estimate. Finally, our results are in agreement with those 
obtained previously by numerical methods [jj, |2| • 

The large K-factors and the significant factorization scale dependence, both related to the opening of new channels 
at NLO, suggest the presence of non negligible NNLO effects. This feature, altogether with the fact that data is 
reasonably described by NLO estimates within uncertainties, obliterate any further effect that might contribute, like 
the resolved component of the photon or those coming from BFKL dynamics in the present experiments. 
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